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Abstract 

Searches for neutral Higgs bosons are performed with the 237 pb~ x of data collected in 
1999 by the ALEPH detector at LEP, for centre-of-mass energies between 191.6 and 
201.6 GeV. These searches apply to Higgs bosons within the context of the Standard 
Model and its minimal super symmetric extension (MSSM) as well as to invisibly decaying 
Higgs bosons. No evidence of a signal is seen. A lower limit on the mass of the Standard 
Model Higgs boson of 107.7 GeV/c 2 at 95% confidence level is set. In the MSSM, lower 
limits of 91.2 and 91.6GeV/c 2 are derived for the masses of the neutral Higgs bosons h 
and A, respectively. For a Higgs boson decaying invisibly and produced with the Standard 
Model cross section, masses below 106.4 GeV/c 2 are excluded. 
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1 Introduction 



Searches for neutral Higgs bosons of the Standard Model, its minimal supersymmetric extension 
(MSSM), and extensions allowing Higgs boson decays into invisible final states were performed 
using data collected by the ALEPH detector at LEP during 1999. The data sample was taken 
at four centre-of-mass energies, 191.6, 195.5, 199.5, and 201.6 GeV at which 28.9, 79.9, 86.3, 
and 41.9 pb _1 of data were collected, respectively. The total data sample corresponding to 
an integrated luminosity of 237 pb -1 was analysed to search for topologies arising from the 
e + e~— > hZ Higgsstrahlung process supplemented by W and Z gauge boson fusion, and from 
the e + e~— > hA associated pair-production process of the MSSM. The production cross section 
of the Higgsstrahlung process in the MSSM is reduced by a factor sin 2 (/3 — a), where tan/3 is 
the ratio of the vacuum expectation values of the two Higgs doublets and a is the mixing angle 
in the CP-even Higgs sector. The hA production cross section is proportional to cos 2 (/3 — a). 
For an invisibly decaying Higgs boson, the observable rate of the Higgsstrahlung process can 
be expressed as £VsM(e + e _ — > hZ), where £ 2 is the product of the branching ratio to invisible 
decays and a model-dependent factor which reduces the cross section with respect to that in 
the Standard Model. 

Searches for neutral Higgs bosons with the ALEPH detector have already been carried out 
up to a centre-of-mass energy of 188.6 GeV [0, 0; no evidence of a signal was found. A lower 
limit at 95% confidence level (CL) was set at 92.9 GeV/c 2 on the Standard Model Higgs boson 
mass. In the MSSM for the benchmark parameter set with maximal stop mixing |J, lower 
limits of 82.5 GeV/c 2 and 82.6 GeV/c 2 were derived on the masses of the h and A Higgs bosons, 
respectively. For an invisibly decaying Higgs boson, a lower limit of 95.4GeV/c 2 was set, for 
a production cross section equal to that in the Standard Model. Similar searches have been 
performed by the other LEP experiments [[|. 

The higher centre-of-mass energies and integrated luminosity in the 1999 data substantially 
increase the experimental sensitivity for the detection of Higgs bosons with respect to previous 
results. Nevertheless, the background processes are the same as those described in Refs. Jl|, |2j. 
The theoretical framework, the event selections, the study of systematic uncertainties, and the 
result extraction are therefore very similar to those previously described. The differences mainly 
consist of reoptimization of event selections, the introduction of new discriminating variables, 
and improvements related to the simulation of signal processes. 

2 ALEPH detector 

The components of the ALEPH detector that are most relevant for the analyses presented here 
are summarized in this section. A more detailed description of the detector can be found in 
Ref. U and its performance in Ref. ||. 

Three coaxial tracking devices are located inside a solenoidal superconducting coil which 
produces an axial magnetic field of 1.5 T. The vertex detector (VDET) J7| consists of two 
cylindrical layers of silicon wafers situated at average radii of 6.3 and 11.0 cm. Charged 
particles with a polar angle in the range |cos#| < 0.88 (0.95) traverse at least two (one) VDET 
layers. The VDET is surrounded by an inner tracking wire chamber (ITC) and by a large 
time projection chamber (TPC), which measures up to 21 three-dimensional coordinates per 
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charged particle between radii of 31 and 180 cm. The tracking achieves a momentum resolution 
<j(pr)/Vr of 6 x 10~ 4 © 0.005, with p^ in GeV/c. The resolution on the three-dimensional 
impact parameter of tracks can be parametrized as (34 + 70/p)(l + 1.6 cos 4 9) /zm, with p in 
GeV/c. 

The electromagnetic calorimeter (ECAL) is also situated inside the coil. It is segmented into 
projective towers of typically 0.9° x 0.9°, which allows electrons and photons to be identified 
within jets. Luminosity calorimeters of similar construction to ECAL are installed between the 
endcaps and the beam pipe and are treated as an extension of the ECAL. A silicon-tungsten 
sampling calorimeter extends the electromagnetic calorimeter coverage down to 34mrad. 

Outside the coil, a hadron calorimeter (HCAL) measures the hadronic energy, acts as a 
filter for the identification of muons, and serves as a return yoke for the magnetic field. The 
outermost detectors are two double layers of muon chambers. 

The measurements of charged particle tracks and of energy deposits in the calorimeters, 
combined with the identification of photons, electrons, and muons, are used to produce a list 
of charged and neutral energy flow particles. Hadron jets, formed by clustering these particles, 
have an energy resolution of cr(E) = (O.QOVE + 0.6) (1 + cos 2 6), where E is the energy in GeV 
and 9 the polar angle of the jet. The resolution on the jet angles is approximately 20mrad in 
both 9 and 0. 

The tagging of b quark jets is accomplished by combining several discriminating 
variables together using neural networks. These variables include an impact-parameter-based 
probability |§, a displaced vertex y 2 an d the transverse momentum of identified leptons 
with respect to the jet axis [|1C ] . Three jet-shape quantities supplement these variables. The 
resulting six variable neural network b-tag is described in Ref. |L1| ; a more specialized four 
variable b-tag used for the e + e~— > hA^ bbbb final state is described in Ref. PJ. 



3 Update of the searches 



3.1 Signal and background simulation 

For each of the four centre-of-mass energies, fully simulated samples of signal and background 
events were generated. The hzha program [[HJ, used to generate the signal events, was also 
used to compute the signal cross sections, the Higgs boson decay branching fractions, and the 
radiative corrections to the Higgs boson masses in the MSSM. The most recent version contains, 
in particular, the latest refinements on these radiative corrections |3|, 13], [L4] and a more complete 
simulation of the he + e~ and hz/ P z/ P final states which now includes the interference between the 



boson fusion and the Higgsstrahlung processes |L5 

The e + e~— > qq(7) background was generated with koralz [TB|, instead of pythia [[T7J, 
to benefit from a more accurate treatment of multiple initial state photon radiation. The 
e + e~— > Weu events, generated with pythia, were compared to events generated with 
and the differences were taken into account in the systematic uncertainties. 
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The other Standard Model four-fermion processes were simulated in a manner similar to that 
described in Refs. |], ||. The sizes of the simulated signal and background samples correspond 
to at least 50 times the collected luminosity used in the analyses. 



2 



3.2 Event selections 



Analyses had already been developed [I], for most final states relevant for the searches for 
neutral Higgs bosons. For the visible e + e - ^ hZ search, these final states are the leptonic final 
state (h£ + £~ where i denotes an electron or a muon), the missing energy final state (huu), 
the four-jet final state (hqq), and the tau final states (hr + r~ and h^ r + r _ , Z— > qq). For the 
e + e~— > hA search, they are the four-b final state (bbbb) and the tau final state (bbr + T~). For 
the invisible e + e~— > hZ search, they are the acoplanar lepton pair final state (h£ + £~) and the 
acoplanar jet pair final state (hqq). As in Ref. fl|], the visible Higgs boson search was conducted 
in two alternative analysis streams, a neural-network-based stream (NN) and a cut-based stream 
(cut), although both streams share the same leptonic and tau-final-state selections. This section 
highlights changes with respect to Refs. |2[ . Most of the hZ and hA selections were optimized 
for Higgs boson masses of 107 and 90GeV/c 2 , respectively, i.e., close to the search sensitivity 
in both cases. 

A reoptimization of the selection criteria and retraining of the neural networks in all channels 
were performed to account for the higher centre-of-mass energies and integrated luminosity. In 
some cases, several new discriminating variables were adopted and refinements to analysis 
techniques implemented, as described below. 

• In the missing energy channel, the reconstructed Higgs boson mass was removed from the 
variables of the single neural network selection (referred to as "A" ) in order to increase the 
selection efficiency for low Higgs boson masses. The structure of the three-neural-network 
selection (referred to as "B") remains unchanged. 

• In the tau final state selection, the overlap with leptonic-final-state events (he + e~ or 
h/i + /i~) was reduced by requiring that the measured invariant mass of the two tau jets be 
smaller than 75 GeV/c 2 when a decay particle of at least one of the two taus is identified 
as an electron or a muon. 

• In the four-jet channels, the following improvements were made. 

i) A four-constraint fit, identical to that used in Ref. [fE|, was adopted in all 



selections, instead of a simple energy rescaling. It allows the jet angles to 
be fitted as well as the jet energies, thus improving the dijet mass resolution, 
particularly near the kinematic threshold. 

In the four-b channel, a new variable was introduced to characterize the 
unbalanced topology of three jets recoiling against one. For each of the possible 
three-jet combinations, the energy-weighted angular dispersion is computed as 




^i£jet i,j,k 

where the indices i, j, and k refer to jets, 9i$ is the angle between the momentum 
direction of the object / contained in jet i, j, k and the vectorial sum u of the 
three jet momenta, and E\ is the energy of object I. The variable A83 is defined 
as the minimum value of this dispersion for all four possible three-jet systems. 
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Figure 1: Distribution, at preselection level, of the minimum value of the energy- weighted angular 
dispersion A#3 in the four-b channel for data (dots with error bars), simulated background (solid 
histogram), and the signal for = ijia = 90GeV/c 2 (dashed histogram). The signal histogram has 
an arbitrary normalization. 



The distribution of this quantity is shown in Fig. [l|. The cut A63 > 50° reduces 
the background from events with gluon splitting e + e~— > bb(g— > bb). 

iii) In the neural network selection, the former 17 variables were supplemented with 
the decay angles relative to the direction of flight of the Z and Higgs boson 
candidates, «i2 and a^. The output of the neural network for the four-jet final 
state is shown in Fig. |2|. 



3.3 Systematic uncertainties 

Systematic uncertainties caused by inaccuracies of the simulation were evaluated as in 
Refs. [1, 2]. Whenever possible, they were extracted from 3.5 pb -1 of data taken at the Z peak 
during the same year. A discrepancy between data and simulation, affecting the acceptance 
for signal and for most of the background processes, was observed in the impact-parameter- 
based b tagging quantities. To correct for this effect, a smearing of the track parameters was 
performed on the simulated events to bring them in better agreement with the data. Half 
of the correction was conservatively taken as a systematic uncertainty for all channels. The 
systematic uncertainties on signal efficiencies are typically 5%, and those on the background 
yields are summarized in Table |1|. 



4 Combination and results 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

NN output 

Figure 2: Distributions of the neural network output used to select four-jet hZ candidates for data 
(dots with error bars), simulated background (solid histograms), and simulated Higgs signal with an 
arbitrary normalization for = 107GeV/c 2 (dashed histogram), after preselection criteria have been 
applied. 



4.1 Selection overlaps 

Events selected by more than one analysis were treated with a well-defined prescription to 
remove overlaps. This procedure prevents any given event from affecting the confidence level 
calculations more than once. The potential overlaps of the searches for different topologies 
arising from the Higgsstrahiung process were suppressed by assigning a precedence order to the 
four selections, and by rejecting any candidate event already selected by a higher precedence 
search. The order of the precedence was chosen to be 1) the missing energy final state; 2) the 
tau final states; 3) the leptonic final state; and 4) the four-jet final state. In the data, no event 
was affected by this procedure. 

The significant overlaps of the hZ and hA searches for tau final states as well as for four-jet 
final states in the cut stream were treated in a manner similar to that used in Ref. |I[ for the 
missing energy channel. This method divides the analyses into three statistically independent 
branches: one branch containing the events selected only by the hZ search; a second branch 
containing the events selected only by the hA search; and a third branch containing the events 
selected simultaneously by both searches. The use of different reconstructed Higgs boson mass 
definitions in the hZ searches (with the Z mass constraint) and the hA searches (with an equal 
mass constraint) introduces an additional difficulty, solved in the following manner. The hZ 
(hA) mass definition is used for the hZ(hA)-oniy branch. For the third branch, the mass 
definition depends on the MSSM parameter set and is chosen such that it maximizes the 
combined expected 95% CL sensitivity. Typically, the hZ mass definition is chosen for small 
tan/3 values, where the hZ cross section dominates, and the hA mass definition is chosen for 
large tan (3 values. 
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Table 1: The numbers of signal and background events expected, and the numbers of candidate events 
observed in the data for each channel. The signal expectation is computed for the hZ process with 
a Higgs boson mass of 107GeV/c 2 , and that for the hA process with a common Higgs boson mass 
of 90GeV/c 2 . The background yield is subdivided into ZZ (including Zee and Zi/P), WW (including 
Wef), and ff (including 77 — > ff). The ff composition is entirely qq for the visible Higgs boson and 
hadronic invisible Higgs boson decay channels, and £ + £~ and 77 — > £ + £~ for the leptonic invisible 
Higgs boson decay channel. Systematic uncertainties are given for each background. 
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1.8±0.2 0.4±0.1 0.1±0.1 


2.3±0.2 





hZ only 


0.3 


0.1 


4.3±0.2 4.6±0.2 0.8±0.2 


9.7±0.3 


13 


hA only 


0.0 


0.0 


0.2±0.0 O.liO.O 0.1±0.1 


0.4±0.1 


2 


invisible h£ + £ 


0.3 




3.4±0.4 3.7±0.5 1.4±0.2 


8.5±0.7 


8 


invisible hqq 


3.3 




31.8±0.8 17.4±3.5 3.1±1.2 


52.3±3.8 


50 



The three-branch subdivision in the NN stream four-jet final state was replaced by a simpler 
and similarly sensitive two-branch treatment: the first branch (4b) containing the events 
selected by the four-b final state search, and the second branch (2b) containing the events 
not selected by the 4b search. 

4.2 Selection results 

For each analysis, the numbers of signal and background events expected, and the observed 
number of events are summarized in Table [I]. For the Standard Model Higgs boson searches, 
a total of 76 events are selected in the data with the NN stream, while 97.8 events are 
expected from Standard Model processes. Figure |3] shows the reconstructed Higgs boson mass 
distributions for the data and the background expectations for the hZ selections. The hA 
selections yield 5 events in the data with the NN stream, while 9.4 events are expected from 
Standard Model processes. In the search for invisibly decaying Higgs bosons, a total of 58 
candidate events are selected with 60.8 events expected from Standard Model processes; the 
distributions of their reconstructed Higgs boson mass are shown in Fig. |. 
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Figure 3: Distributions of the reconstructed Higgs boson mass for the data (dots with error bars) 
selected in the hZ searches by (a) the NN stream, and (b) the cut stream. The histograms show the 
Standard Model background expectation. 



Fewer events were selected in the data than were expected. This deficit is mostly due to 
the four-jet analysis of the NN stream where 46.4 events were expected, but only 30 events 
were selected in the data. In the four-jet cut analysis 27 events are selected in data with 
36.3 expected. As can be seen in Fig. j|, the deficit of events is mostly near the mass region 
dominated by e + e - -^ ZZ events. 

At the preselection level, the background expectation is dominated by four-jet 
e + e~— > W + W~ events which are kinematically similar to signal events, and no significant 
discrepancies were found between data and Monte Carlo simulation for all relevant kinematic 
variables. Figure [| shows the distribution of the neural network output distribution at the 
preselection level. Overall, the agreement between the simulated expectation and the data 
is very good with 2572 events expected and 2523 seen in the data. The deficit is, however, 
apparent in Fig. ^| at high values of the neural network output, which is dominated by b tagged 
e + e~— > TIL events. 

As mentioned in Section 3~7 3", a slight disagreement was found and corrected for in the 
b tagging distributions. The corrections were determined from Z peak data taken during 
the same year. These were verified with e + e~— > qqy radiative return events and semileptonic 
e + e~— > W + W~^ qqiu events collected at high centre-of-mass energies. The qcfy events provide 
a known bb content while the semileptonic W + W~ events provide a control sample with 
negligible b-jet content. The impact of these b-tagging corrections is propagated to the 
simulation, but the effect is too small to account for the observed deficit of events. The 
conclusion of these studies is that the cause of the deficit is most likely due to a statistical 
fluctuation. 
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Figure 4: Distribution of the reconstructed Higgs boson mass Mh for the data (dots with error 
bars) and the expected background (histogram) in the invisibly decaying Higgs boson searches for the 
acoplanar lepton and the acoplanar jet final states. 



4.3 Exclusion limits 

The 176.2 pb -1 of data collected in 1998 with a centre-of-mass energy of 188.6 GeV were 
combined with the 1999 data to determine the final results. As was done in Refs. [p], [2|, 
confidence levels on the signal hypothesis were drawn using the Signal Estimator method pD 



rather than with the method used by the other LEP experiments ||21|| , thus increasing the 
expected mass limit by ~ 0.5 GeV/c 2 . Expected confidence levels were calculated as the median 
confidence level in the absence of signal. 

The tau final state analyses and the cut analyses for the four-jet, four-b, and missing energy 
topologies use only the reconstructed Higgs boson mass as a discriminant in the confidence level 
calculations. The other analyses use an additional discriminant variable: the neural network 
output in the case of the 2b branch of the NN stream four-jet analysis and a b-tag sensitive 
variable in the 4b branch as well as in the h£ + £~ final state analysis. 

The expected and observed confidence levels are shown in Fig. [| for the Standard Model 
Higgs boson searches. The low confidence levels for the background hypothesis at small Higgs 



boson masses are a reflection of the observed deficit of events discussed in Section J4.2I. The lower 



limits at 95% confidence level on the mass of the Standard Model Higgs boson with systematic 



uncertainties taken into account p2| are summarized in Table Since the NN stream gives an 
expected Higgs boson mass lower limit higher by about 0.5 GeV/c 2 , it is used to set the final 
results. 

Figure H shows the hZ, the hA, and their combined Higgs boson mass exclusions as a function 
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Table 2: Summary of the 95% CL lower limits on Higgs boson masses in the Standard Model and 
in the MSSM, and their expected values in absence of signal. The tan f3 excluded ranges in the m™ ax 
scenario are also given for a top quark mass of 175GeV/c 2 . 





95% CL lower limits (GeV/c 2 ) 




NN stream 


Cut stream 




Median Expected 


Observed 


Median Expected 


Observed 


SM m h 


107.8 


107.7 


107.3 


105.2 


MSSM m h 


88.7 


91.2 


88.3 


90.7 


MSSM m A 


89.1 


91.6 


88.7 


91.1 


tan (3 


[0.8,1-9] 


[0.8,1-9] 


[0.8,1.9] 


[0.8, 1.7] 



of sm 2 ((3 — a). These combined exclusions are interpreted within two benchmark parameter 
sets, one which maximizes the radiative corrections to the lighter CP-even Higgs boson mass 
as a function of tan (3 (referred to as the m™ 5 ™ scenario) JT^j and one in which no mixing in the 
scalar top sector is assumed. The excluded domains in the [m^tan/?] plane are shown in Fig. [7| 
as obtained using the NN stream. The results in the m™ ax scenario are summarized in Table 0. 
The tan /3 exclusion range is sensitive to the top quark mass assumption which was taken to be 
175 GeV/c 2 . For a top quark mass of 180 GeV/c 2 , the tan/3 exclusion range reduces to [1.0,1.5]. 

For the invisibly decaying Higgs boson search, no discriminant variable is used in the 
confidence level calculations, as was done in the previous search at 188.6 GeV @, because 
the sliding analysis technique effectively acts as if the reconstructed mass were used as a 
discriminant variable. The results are interpreted as an exclusion domain in the (mh,( 2 ) plane, 
presented in Fig. || For £ 2 = 1, invisibly decaying Higgs bosons with masses below 106.4 GeV/ c 2 
are excluded at 95% CL with an expected limit of 105.6 GeV/c 2 . 

5 Conclusions 

Searches for neutral Higgs bosons in e + e~ collisions at centre-of-mass energies from 191.6 to 
201.6 GeV have been carried out with the ALEPH detector using an integrated luminosity 
of 237 pb~ x . The major event topologies expected from the hZ process for visibly as well as 
invisibly decaying Higgs bosons, the weak boson fusion process, and the hA process have been 
analysed. The data were combined with the 176 pb -1 collected in 1998 at a centre-of-mass 
energy of 188.6 GeV. 

A 95% CL lower limit on the Standard Model Higgs boson mass of 107.7 GeV/c 2 is obtained 
with an expected sensitivity of 107.8 GeV/c 2 . Lower limits of 91.2 and 91.6 GeV/c 2 are obtained 
at 95% CL for the masses of the h and A neutral Higgs bosons of the MSSM. These limits, 
determined in the m™ ax scenario, are valid for any tan/3 > 0.5. For a top quark mass of 
175GeV/c 2 , the tan/3 range between 0.8 and 1.9 is excluded with 95% confidence. 

An invisibly decaying Higgs boson with a production cross section equal to that in the 
Standard Model is excluded at 95% CL for masses below 106.4 GeV/c 2 . 
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Figure 5: Observed (solid) and expected (dashed) CL curves for (a) the background hypothesis and 
(b) the signal hypothesis for the NN stream, as functions of the hypothesized Standard Model Higgs 
boson mass m^. For the signal hypothesis (b), the intersections of the horizontal line at 5% with the 
curves define the observed and expected 95% CL lower limits on the Standard Model Higgs boson 
mass. 




Figure 6: The expected (dashed) and observed (solid) excluded regions at 95% CL in the plane 
[mh,sin 2 (/3 — a)] for (a) the NN stream and (b) the cut stream. 
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Figure 7: The expected (dashed) and observed (solid) excluded regions at 95% CL for the NN stream 
in the [mh, tan/?] plane of the MSSM for (a) the m™ ax parameter set and (b) the no mixing scenario. 
The densely hatched regions are not allowed theoretically. 
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Figure 8: Observed (solid) and expected (dashed) regions excluded at the 95% CL in the (m^,^ 2 
plane for an invisibly decaying Higgs boson. 
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